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After a brief discussion of the promise and limitations of the lattice technique, I review lattice QCD results 
for several quantities of phenomenological interest. These are: matrix elements for heavy-to-light meson decays, 
leptonic decay constants fg and / jj , the parameters Bb and Bk for neutral B and K meson mixing respectively, 
the strong coupling constant, light quark masses and the lightest scalar glueball mass. 



1 Introduction 

In this review I concentrate on lattice QCD results 
for a selection of quantities directly relevant for phe- 
nomenology, some of which were not otherwise re- 
ported at this conference. I apologise to presenters 
of lattice studies which I do not have space to report 
on here and refer the interested reader directly to the 
parallel session reports]] 

Lattice QCD is an important tool for the non- 
perturbative evaluation of strong interaction effects, 
but a wary consumer should keep in mind the main 
sources of error, which I will briefly mention below, 
when using lattice results. One beauty of the lattice 
approach is that these errors can be systematically 
investigated and reduced. 

The presentation begins with form factors for 
B°^ir+l-vi, B°^p + l-Di and B^K*~f. The semilep- 
tonic decays are important for determining \V u b\] the 
radiative decay for Vt s and as a window on new 
physics. Next follow results for the B meson decay 
constant fs and mixing parameter Bb, crucial for 
constraining the Standard Model unitarity triangle. 

1 then turn to the kaon mixing parameter Bk, strong 
coupling a s and the light quark masses, concluding 
with results for the lightest scalar glueball. 

2 Lattice Calculations 

The standard lattice approach to QCD uses a discre- 
tised finite volume region of Euclidean space-time on 
which the quantum field theory path integral becomes 
a well-defined multi-dimensional integral evaluated by 
Monte Carlo methods. 

Current simulations use lattice spacings a in the 
range 0.05-0.2 fm, sufficient to cover a hadron of size 
0.5 fm or greater by a few lattice points. These a's cor- 
respond to energy scales of 1-4 GeV putting the lattice 
ultraviolet cutoff above the scale of low energy QCD 
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dynamics. Continuum results should be obtained in 
the limit a— >0: this continuum extrapolation is now 
feasible for many quantities. Lattices should be large 
enough that hadrons will comfortably fit on them. 
Current spatial sizes are of order 2 fm or so, though 
some practitioners advocate at least 2.5fmlil 

Because the QCD action is quadratic in the quark 
fields, matrix elements of quarks can be evaluated us- 
ing quark propagators in the gluon background to- 
gether with the gluon-field-dependent determinant of 
the fermion operator. The determinant is extremely 
demanding to calculate so it is often set to its average 
value in the gauge field background — the quenched 
approximation which corresponds to neglecting in- 
ternal quark loops. In practice, much of the effect 
of internal loops is to change the running of the cou- 
pling constant and this can be compensated by chang- 
ing the value of /?, the bare coupling which is input. 
Quenching is one of the systematic effects causing dis- 
agreement between lattice spacings determined from 
different physical quantities. 

Calculating quark propagators means inverting 
the fermion operator. This is very slow for realistic 
light quarks, so a range of masses around the strange 
mass is simulated and then a 'chiral extrapolation' 
made to realistic values. For heavy quark masses m, 
the inversion is fast but for large ma discretisation er- 
rors are large. Hence, B physics results using relativis- 
tic fermion actions typically involve an extrapolation 
from results at masses close to the charm scale. An al- 
ternative is to use static or nonrelativistic (NRQCD) 
lattice actions. In the static case the quark mass is 
treated analytically, outside the simulation, and re- 
sults are obtained in a systematic expansion around 
the infinite mass limit. NRQCD actions work well for 
quarks around the b mass and above, but begin to fail 
as one approaches the c mass. Actions are also used 
which interpolate between the static and relativistic 
extremesa 

Lattice calculations provide matrix elements of 
bare operators defined with the lattice regularisation, 
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but results are needed for matrix elements in a contin- 
uum scheme, like MS. A continuum operator typically 
matches onto a set of lattice operators, 

Of>»\n) = Z ij {n,a)OY\a), 

where the renormalisation constants Zij can be cal- 
culated perturbatively. Lattice perturbation theory 
using the bare coupling t/ 2 , = 6//3 is notoriously badly- 
behaved, however, u Abandoning g 2 . in favour of a 
continuum-like definition leads to better behaviour, 
but nonperturbative methods are b^ing developed and 
used (some are mentioned below) cl This looks to be 
the way of the future. 

Lattice gluon actions have discretisation errors 
beginning at 0(a 2 ). The simplest, Wilson, fermion 
action, a latticised relativistic Dirac action, has errors 
beginning at O(a). One can simply accept the O(a) 
errors and let the a— >0 extrapolation remove them. 
However, removing or reducing the 0(a) errors allows 
simulations at coarser lattice spacings and makes the 
continuum extrapolation less severe. This 'improve- 
ment' of lattice actions is currently.-|a very active sub- 
ject. The Sheikholeslami-Wohlert □ (SW or 'clover') 
action removes O(a) errors at tree level, so the first 
corrections are 0(a s a). The SW action contains one 
extra term compared to the Wilson action. By tuning 
the coefficient of this term one can reduce or even re- 
move all 0(a) errors in physical quantities (one has to 
improve operators used in matrix element calculations 
also). Partial removal using 'tadpole-improvement' 13 
is being widely applied, but-an ambitious program by 
the ALPHA collaboration ErtH aims for complete re- 
moval. 

3 Semileptonic and Radiative Heavy-to-light 

Decays 

Lattice form factor calculations are crucial here: the 
overall normalisation at the zero recoil point uj = 
v-v' — 1 is not fixed by heavy quark symmetry as 
it is for a heavy-to-heavy transition. Here v and v' 
are the four- velocities of the B meson of mass M and 
the meson of mass m it decays into, respectively. The 
squared momentum transfer to the leptons or photon 



is q z = M z 



2Mmu. 



Relativistic quark calculations use heavy quark 
masses around the charm mass. The initial heavy me- 
son is given or 1 lattice units of three-momentum, 
while the light final meson can generally be given up 
to two lattice units of spatial momentum, allowing q 2 
to be varied from g^ ax = (M— m) 2 (where uj = 1) 
down to q 2 < at the D scale. Heavy quark symme- 
try determines the M dependence of the form factors 
at fixed u>, enabling an extrapolation to the B scale. 



g 2 /GeV 




M/GeV 



Figure 1: q 2 range for heavy-to-light decays as a function of the 
decaying heavy meson mass. Lines of constant ui are shown. 
The light final state mass is taken to be 850 MeV, typical of 
lattice pseudoscalar or vector meson masses before chiral ex- 
trapolation. 



Fig. [I] shows that scaling in M at fixed u> sweeps all 
the measured points to a region near g^ ax for B de- 
cays. The problem is then to extrapolate back down 
to q 2 — 0. This is particularly acute for the radiative 
decay B^K*^ where only the form factors at q 2 = 
contribute. Even if a static, non-relativistic or other 
modified action is used for the heavy quark, the re- 
striction on usable three-momenta on current lattices, 
caused by momentum-dependent errors and increas- 
ing statistical uncertainty, ensures., that form factor 
values are obtained only near g^xul 

Some assistance is provided by ensuring that any 
model q 2 dependences respect heavy quark symme- 
try and constraints relating form factors at q 2 = 0. 
For example, / + and f° for B — ► it are related at 
q 2 = and consistency can be achieved by fitting f + 
to a dipole [pole] form and /° to a pole [constant] 
form. Heavy quark symmetry and light flavour SU (3) 
symmetry relate the B°^p + l~i>i and B—^K*^ form 
factors. Models can further relate these to form fac- 
tors for B°^tt + 1~Di. An overall fit might then be 
used. However, it is clearly desirable to avoid models 
entirely. 

3. 1 Semileptonic B — ► 7r 



Lattice results Bffla for the B° — >tt + 1~Di form fac- 
tors, p£|t-| and /°, are plotted in Fig. || (further re- 
sultsli§t£l at q 2 = only are not displayed). 

For massless leptons, the decay rate is determined 
by / + alone. However, the constraint / + (0) = /°(0) 
(with suitable conventions) makes lattice measure- 
ments of both form factors useful. One procedure uses 
dispersive constraints, obtained by combining disper- 
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Figure 2:r-fiattice measuremapta of the form factors for B°— >7r+i ui and B°—>p + l i/;.rRoints are from ELcE (=), APE00 (=), 
UKQCD £3 (•) and FNALEE3 (x, preliminary). Systematic errors have been added fed to the UKQCD and FNAL points. The 

vertical dashed lines show (jfmax- 



sion relations with unitarity, analyticity, crossing sym- 



,to obtain model inde- 
has shown how to in- 



metry and perturbative QC 
pendent bounds. Lclloucht 
corporate the ,/ + (0) = /°(0) condition together with 
imperfectly known values of the form factors, typi- 
cal of lattice results with errors, to obtain families of 
bounds with varying confidence levels. A set of such 
bounds is shown in Fig. || together with the UKQCDLL3 
lattice measurements used to obtain them. In the fig- 
ure, f° and / + are plotted back-to-back, showing the 
effect of imposing the constraint at q 2 =0. 

In Table |l| the bounds have been used to give 
ranges of values for the decay rate T(B ^tt + 1~Di) in 
units of |V^b| 2 ps _1 together with values for the form 



mental result for the decay rate, one can extract \V u b\ 
with about 35% theoretical error. Although not very 
precise, this result relies only on lattice calculations of 
matrix elements and heavy quark symmetry, together 
with perturbative QCD and analyticity properties in 
applying the dispersive constraints. There is no model 
dependence. Improved lattice results can be used as 
input for the bounds once they become available. 

Also shown in Table |l|, for comparison, are val- 
ues obtained from lattice calculations where assumed 
q 2 dependences have been imposed. For ELC □ and 
APE(L3 one value of /+ has been used, at the given 
value of q 2 , fitted to a singk pole form with pole mass 
m p . The UKQCD result E3 is obtained from a corn- 



factor / + at q 2 = 0. When combined with the experi- bined dipole/pole fit to all measured f + /f° points. 
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Table 1: B°— >7r+Z~ «7-iresults from dispersive constraints ap- 
plied to lattice datafj together with results obtained using 
ansatze for the form factor /+. The decay rates are values 

for rogo^TT+rpo/lv^pps- 1 





Rate 


/ + (0) 


Dispersive _ . 
Constraint li-3 


2.4-28 
3.6-17 
4.8-10 


-0.26-0.92 95% CL 
0.00-0.68 70% CL 
0.18-0.49 30% CL 


ELci 9 ± 6 

q 2 -18 GeV 2 , pole fit, m F 


0.10-0.49 
=5.29(1) GeV 


APEEl 8 ± 4 0.23-0.43 
q 2 ~20.4GeV 2 , pole fit, m p =5.32(l) GeV 


UKQCD Ht5 7±1 
dipole/pole fit to f+/f° 


0.21-0.27 




7 2 /GeV 2 



Figure 3: Bound* on f + and f° for B°— >7r + / — P; fromi-disper- 
sive constraints^ The data points are from UKQCDt3 with 
added systematic errors. The pairs of fine curves are, outer- 
most to innermost, 95%, 70% and 38% bounds. The upp«r and 
lower shaded curves are light-cone 113 and three-point E9 sum 
rule results respectively. 



Note that the UKQCD points have statistical errors 
only and have not been chirally extrapolated — they 
correspond to a pion mass of around-SOO MeV (a sim- 
ilar caveat applies for the F.NALtrt 3 ] results for /+ 
and /°). The results given 113 have used these val- 
ues as though they applied to the physical pion. Lr± 
obtaining bounds based on these points, Lellouchli-3 
added a conservatively estimated systematic error in- 
cluding terms to account for the chiral extrapolation 
(this error has been added to the UKQCD and FNAL 
points plotted in Fig. ||). 

Chiral extrapolations are severe for the B — ► 7r 
matrix elements. As the pion mass approaches its 



physical value the £?*-pole and the beginning of the 
Btt continuum approach <7^ ax from above and the 
form factors may vary rapidly with the pion mass near 
9max- This is not a problem for B — > p so lattice cal- 
culations of semileptonic B — > p decay are currently 
most reliable and I now turn to them. 

3.2 Semileptonic B — ► p 

To avoid models for the q 2 dependence of the form 
factors, we use the lattice results directly. The lattice 
can give the differential decay rate dT/dq 2 , or the par- 
tially integrated rate, in a q 2 region near q 2 nax up to. 
the unknown factor |K,f,| 2 . For example, UKQCD Eil 
parametrised the differential decay rate near q^ ax by, 



-2 = const|Kl,|VA 1/2 a 2 (l + b[q 2 ~q 2 max }) , (1) 

where A is the usual phase space factor and a and b are 
constants. The constant a plays the role of the Isgur- 
Wise function evaluated at u) = 1 for heavy-to-heavy 
transitions, but in this case there is no symmetry to 
determine its value. 

For massless leptons, the differential decay rate 
depends on the V, A\ and A 2 form factors, but A\ 
is the dominant contribution near q^ aK and is also 
the best measured on the lattice, as shown in Fig. ||. 
UKQCD perfcppcd a fit to the parametrisation in 
Eq. (0) to findP 

a = 4,6t° i ± °- 6 GeV - 

Discounting experimental errors, this result will al- 
low determination of \V u b\ with a theoretical un- 
certainty of 10% statistical and 12% systematic. 
CLEO are beginning to extract the differential decay 
distributions Ej 

The UKQCD results for V, A\ and A 2 agree very 
well with a light cpae sum rule (LCSR) calculation 
of Ball and Braunc3 More interestingly, LCSR cal- 
culations predict that all the form factors for heavy- 
to-light decays hay£_the following heavy mass depen- 
dence at q 2 = 0:lijC3 



/(0) 



1 



M 3 / 2 



(a + a 1 /M + a 2 /M 2 + ■■■). (2) 



The leading M dependence comprises yM from the 
heavy state normalisation together with the behaviour 
of the leading twist-light cone wavefunction. For A\ (0) 
the LCSR result isN A x (0) ~ 0.26. UKQCD fitted the 
heavy mass-dependence of A\ (0) to the form in Eq. (||) 
and found I 2 -!! 



^i(0) 



0.18 ±0.02 

n 99 +0.04 
u - zz -0.03 



ao and a\ ^ 
qq, ai and a 2 ^ 
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Table 2: Lattice results for B^K*^. 



On the lattice one calculates a dimensionless quantity 
Zl-, from which fp is obtained via 



T(0) 









M-V2 


r 2(g max ) 


BHSE 




0.10(3) 
0.09(1) 
0.09(1) 

o.i5(D 




0.33(7) 

0.23(3) 
0.27(f) 
0.30(3) 


LANL. 

APEE 

ukqjc 

BHSE 




0.24(1) 
0.23(3) 
0.26(f) 



Pole fits for A\ have leading M 3 / 2 behaviour at q 2 = 
0, so we can also compare with other lattice results: 

f 0.22 ±0.05 ELcl, 
Ai{0) = \ 0.24 ±0.12 APEt3 

[ 0.27±I UKQCDEl 

3.3 Rare radiative B — > iiT*7 

This decay_was discussed in some detail by A. Soni at 
Lattice 95 El so my comments will be brief. Table ^| 
summarises available lattice results, all from quenched 
simulations, for the matrix element 

(K*(k,r,)\sa^q v b R \B(p)) (3) 

which is parameterised by three form factors, Ti, i = 
1,2,3. For the decay rate, the related values Tj (0) and 
T 2 (0) are needed. Suitably defined, they are equal, so 
the table lists a single value T(0), together with the 
directly measured ?2(<7max)- The results are classified 
by the leading M dependence of the form factor at 
q 2 = which is governed by the model used for the q 2 
dependence. Dipole/pole forms for T1/T2 give M~ 3 / 2 
behaviour and pole/constant forms give M -1 / 2 . The 
results agree when the same assumptions are made. 
All groups find that T2 has much less q 2 dependence 
than Ti, but the overall forms cannot be decided, so 
a phcnomcnological prediction is elusive. 

Additional long distance contributions may not 
be negligible so the matrix elfimeirt of Eq. (|3j) may 
not give the true decay ratec^lEl'EJ Once the q 2 de- 
pendence of the form factors is known, lattice calcu- 
lations of the ratio R K * = T(B^K*j)/T{b -> sj) 
can be compared to the experimental result to test 
for long distance effects. 

4 Leptonic Decay Constants Jb and fjj 

The leptonic decay constant of a pseudoscalar meson 
P is defined by the axial current matrix element 

(0\A fl (0)\P{p))=if P p fi . 



f P ^M^~2^Z^Z L a- 3 ' 2 

where Z Ten is the renormalisation constant required to 
match to the continuum and a is the lattice spacing. 

Two collaborations have new values, shown in ta- 
ble [|, for Jb and fp> obtained from continuum extrap- 
olations p£ quenched results at several lattice spacings. 
JLQCD E3 study different prescriptions for reducing 
the 0(ma) discretisation errors associated with heavy 
quarks and aim to show that all results converge in the 
continuum limit. Their error is statistical combined 
with an error from the spread over prescriptions. 

MILC El simulate a range of heavy quark masses 
plus a static (infinite mass) quark, giving meson 
masses straddling mo and allowing an interpolation 
to ma. Their (preliminary) results give errors from: 
(i) statistics, (ii) systematics within the quenched 
approximation and (iii) unquenching — they have dy- 
namical fermion results but do not yet perform a con- 
tinuum extrapolation with them. Their results sug- 
gest that unquenching will raise the value of t.Vie de- 
cay constants. This agrees with an estimatelp'Ea us- 
ing the difference between chiral loop contributions in 
quenched and unquenched QCD, that fs s in full QCD 
is increased by 20% over its quenched value. Calcula- 
tions of Jb in the static limit extrapolated from nega- 
tive numbers of flavours c3 also suggest an increase of 
about 20%. 

The SGO collaborationcilo is calculating fs us- 
ing a lattice NRQCD action (and heavy-light axial 
current) corrected to 0(1/M) where M is the heavy 
quark mass. However, the renormalisation constants 
required to match onto full QCD are not yet available, 
so I will not quote results for the decay constant. 

The last row of table || summarises £□ global re- 
sults for B and D meson decay constants following 
the compilation by G. Martinelli.tjI^'LD 

5 B—B Mixing Parameter 

Recent lattice calculations of the B meson mixing pa- 
rameter Bb (a*) j defined analogously to the kaon mix- 
ing parameter Bk in Eq. ^ below, have been made 
both-ihr- xelativistic quarks and_in_the static 

limit E-a'CJ'E3 The latest static results c§c£l use new 
calculations of the full-theory/static-theory matching 
incorporating previously omitted contributions! 47 !! 48 ! 
Since Bb(h) is scale-dependent it is conventional to 
quote a renormalisation group invariant (RGI) quan- 
tity, Bb- At one loop B B = a s (ny 2/ ^.E^(p), but 
the two-loop relation is commonly usedEJ'E3 
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Table 3: VaU*«s and ratios of leptonic decay constants of pseudoscalar mesons from recent continuum extrapolations by MILCE3 
and JLQCDEj. Statistical and systematic errors have been combined in quadrature. The MILC rea ulta j pa hwle a systematic error 
for quenching. The bottom row contains summary values based on global data EjcjjKr E3 







/s/MeV 


WMeV 


fs a 1 Ib 


f D / MeV 


f D JMeV 


Id 3 1 Id 


milcIe 

JLQCD 




166(33) 


181(41) 
197(a?) 


1.10(10) 


196(18) 
202(H) 


211(28) 
216(?i) 


1.09(7) 


Summary 


175(25) 


200(25) 


1.15(5) 


205(15) 


235(15) 


1.15(5) 



Calculations with relativistic quarks show no ob- 
vious lattice spacing dependence and cluster around 
a value of B B = 1-3(1) [with a weighted average 
of 1.31(3)_janoted in Warsaw] for the one-loop B 
parameter^Jo or Bb = 1-4(1) using the two-loop 
formulaEj Allowing for uncertainty in the continuum 
extrapolation, I quote the two-loop Bb = 1-4(1) as 
the quenched lattice result. 

j-Ihe static, results differ for calculations using WiL 
sonc3 or SWe5e3 light quarks. The Kentucky grouplla 
(Wilson) find B B = 1.40(6) for the one-loop RGI 
quantity, rising to 1^*0(6) for the two-loop case. The 
APE two-loop results (SW) is B B = 1.08(6)(8). The 
dominant uncertainty is from higher order terms in 
the perturbative matching to continuum QCD. Non- 
perturbative renormaJisation will be crucial to reduce 
systematic errorsEJO 

The relevant quantity for B-B mixing is £ B = 
f B B B - Taking the two-loop relativistic quark result 
for B B with f B = 175(25) MeV from table | gives 
£b = 207(30) MeV as the current lattice estimateEB 
This quantity can also be extracted directly from the 
AB=2 matrix element, M = 8f B B B m B /3. To avoid 
uncertainties from setting the scale, it is convenient 
to determine the ratio c3 



Tsd 



pB a B B s m B3 
f B B B m 2 B 



For relativistic quarks a direct extraction gives r sc i = 
1.54(13)(32). Previously, r s d has been evaluated from 
separate results for the decay constant and B param- 
eter ratios. Combining £a s / fs = 1.15(5) from table || 
with BbJBb = 1.01(4)Eland m Bs /m B p. 1.017 gives 
r sd = 1.38(13). In the static case, APEB find r sd = 
1.35(5) by the direct method, or r s d = 1-43(7) from 
combining f B Jf B = 1.17(3) and B B JB B = 1.01(1) 
measured on the same gauge configurations. The re- 
sults of the two methods are quite consistent, but fu- 
ture calculations should improve on the precision of 
r s d obtained directly. 

The above results are from quenched calculations. 
Unqucnching is expected to increase f Bg / f B by about 



10% (0.16 from a chiral loop estimate BE). For 
B Bb /B b , numerical evidence suggests a small increase 
on two-flavour dynamical configurations!^! but the chi- 
ral loop estimateEStJ is for a decrease of —0.04 in the 
ratio. 

6 Kaon S-Parameter Bk 

The parameter Bk is defined by 



8{K \s-ff,Ld\0) {Q\s^Ld\K)/^ 



(4) 



where L = 1— 75. It is a scale dependent quantity for 
which lattice results are most often quoted after trans- 
lation to the value in MS using naive dimensional reg- 
ularisation (NDR) at a scale fj, = 2 GeV. I will follow 
this practice while discussing the lattice results and 
convert at the end to the renormalisation group invari- 
ant parameter Bk normally used in phenomenology. 
At next-to-leading order, Bk is given by, 



B 



K 



Bk{v) 



"s(m)^( 



32. / , a s {fi) (/?i7o -/3o7i) 



4tt 



where /?o.i and 70,1 are the first two coefficients 
of the beta-ipiunction and anomalous dimension, 
respectivelyeaH 

Systematic errors in Bk calculations are be- 
ing carefully explored. Calculations using staggered 
fcrmions (an alternative formulation of relativistic lat- 
tice fermions) are statistically more precise, but Wil- 
son fcrmion results are rapidly improving. Here I sum- 
marise the current situation. For mate details see the 
report by S. Sharpe from Lattice 96E3 

6. 1 Staggered Bk 

Discretisation errors for JLf using staggered fermions 
are known to be 0(a 2 )E3 The first calculation E3 
performed with a range of lattice spacings a there- 
fore used a quadratic extrapolation in a and found 
B K BR (2GeV) = 0.616(20)(27) for the quenched re- 
sult. The data itself, however, could not distinguish 
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Figure 4: JLQCD resultsEl for Bg DR (2GeV) as a function of 
m p a, with a quadratic fit to four leftmost data points (solid) 
and a linear fit to five leftmost data points (dashed). "Non- 
inv" and "Inv" refer to two possible discretisations of the lattice 
operators. 

linear and quadratic a dependence. New results from 
the JLQCD collaboration are shown in Fig. [|. 
Their data fits better to a linear than a quadratic 
dependence on a. The .axgiiaieiit for 0(a 2 ) correc- 



tions has been checked 
quoted for a quadratic 



3' 5 



at 



i3 however, so results are 
B Future calculations at 



smaller a should confirm the leading a 2 dependence. 
A continuum-extrapolated quenched result has also 
been given by a group from The new results 

are (n/ is the number of dynamical flavours): 



D 



NDR 
K 



(2 GeV) a^o 



_ J 0.587(7)(17) 
0.573(15) 



JLQCD 

osu 



(5) 



I will take these as the best quenched lattice estimates 
of B K . 

Two issues need addressing to relate the results 
in Eq. (||) to Bk for full QCD. One is the inclusion of 
dynamical quarks. The second is to allow m s ^ md, 
since the calculations above have a kaon composed 
from degenerate quarks. 

Important progress inj-unquenching Bk has been 
made by the OSU grouppa who find a statistically 
significant increase in Bk (earlier studies suggested a 
small decreased), 



dNDR 
a K 



(2 GeV)' 



fl£ DR (2GeV) 



n,=0 



1.05(2). 



(6) 



They calculated with u/ = 0, 2, 4 for fixed lattice spac- 
ing a — 2 GeV. Sharped gives a more conservative 
estimate for the ratio as 1.05(15) after allowing for the 
extrapolation to a=0. 



To calculate Bk with physical mass non- 
degenerate quarks requires fermions with very small 
masses and will therefore be difficult. However, chi- 
ral perturbation theory fixes the quark mass depen- 
dence of Bk, providing. 2l. way to determine the non- 
degeneracy correctipnl£3'E3 The calculation has yet to 
be done, so Sharped estimates, 



B 



nondegen 
K 



1.05(5) B™ 



degen 



(7) 



Combining the staggered fcrmion results in 
Eq. (^) with the unquenching and nondegeneracy cor- 
rections in Eqs. (||) and (^) leads to the final estimate: 



dNDR 



(2 GeV) = 0.64(2)(10). 



(8) 



The first error is that in the quenched value. The 
second is the larger 15% unquenching error combined 
in quadrature with the 5% error for nondegeneracy. 
Since Eq. (||) incorporates estimates of systematic ef- 
fects in the central value, an alternative statement is 
B K DR (2GeV) = 0.58(2)(9), where the central value 
is the quenched result, noting that unquenching and 
nondegeneracy can raise the value by 10%. Convert- 
ing the result in Eq. (||) to B K using a s (2GeV) = 0.3 
and three flavours gives || 

Bk — 0.87(3)(14). 

6.2 Wilson B K 

Calculations of Bk using Wilson and SW fermions 
have to deal with the explicit breaking of chiral sym- 
metry by the fermion action. This means that the 
continuum operator of interest mixes with four other 
dimension six lattice operators: 



The constants Z and Zi have to be determined so that 
lattice and continuum matrix elements agree to O(a). 
A general four-fermion operator has the chiral expan- 
sion 



(K\0\K) = a 
Pk-Pk{j- 



\-f3m 2 K 
52m 2 K 



A- 



- 8\mj 



For the continuum operator used to determine Bk, 
chiral symmetry demands that a — (3 = 8\ =0. The 
vanishing of these momentum-independent lattice ar- 
tifacts can be used to test that the Zi are correct. 



a This-*ffers from 0.90 
JLQCD E3 and new OSU 
unquenching and degeneracy error 



quoted in Warsaw. It uses updatpdi 
results with a more conservative E3 
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Various methods have been applied to determine 
the Zi and Z . One loop perturbation theory gives in- 
correct chiral behaviour. One can adiu&Lthe z% by 
hand to restore the correct behaviour or calcu- 
late with ya.rving momenta to isolate and discard the 
artifactsEjEj A better procedure is to calculate Z and 
the Zi nonperturbatively. The Rome group demand 
that quark matrix clpmiaits satisfy contimium nor- 
malisation conditionsproQ while JLQCD E3 impose 
chiral Ward identities on quark matrix elements to 
determine the z%, together with a continuum normali- 
sation step to fix Z . Both groups find that the matrix 
element of O cont is obtained with the correct chiral 
behaviour. The errors are currently larger than for 
staggered fermions, principally because of the need to 
calculate several constants Zi. 

JLQCD find that the continuum extrapolation 
is best done for a quantity differing from Bk by 
lattice artifacts which vanish as a— +0. Extrapolat- 
ing from results at three different a values they find 
££ DR (2GeV)™£^° = 0.59(8), in excellent agreement 
with the staggered results in Eq. (||). 

7 Strong Coupling a s 

Determinations of a s from Lattice QCD are done in 
three steps: (i) define and measure some a™**, (ii) 
determine the lattice spacing a, to set the scale at 
which a' att takes its measured value, and (iii) convert 
to tt^jg. Step (i) is necessary because the bare lattice 
coupling, determined from the simulation parameter 
(3 = 6/g 2 , is rather small and has a badly-behaved 
perturbation theory. One must instead use a phys- 
ical definition for a|, att . Steps (ii) and (iii) are the 
major source of uncertainty. Here I will give an up- 
date andj-reJer the reader to recent reviews for more 
detailsE3i£3'E3 New results are available 123 from the 
NRQCD collaboration and a Fermilab-SCRI group, 
both using quarkonium level splittings to fix the lat- 
tice spacing. Other lattice methods of determining a s 
are being developed. These include studying the three 
gluon vertex EJ and a program by the ALPHA col- 
laboration using the Schrodinger Functional method, 
outlined by S. Sint at this conference! 




5 10 

Q = 3.4/a (GeV) 

Figure 5: Recent lattice determinations of l/«p as p-function 
of the scale Q = 3.4/a at which they are measured Lj. Boxes 
denote quenched (n^=0) results with the scale set from char- 
monium and circles are quenched results with the scale from T. 
Two-flavour (ny=2) results from charmonium and T are de- 
noted by pluses and crosses respectively. Results extrapolated 
to i/=3 are given by fancy boxes. The upper and lower solid 
curves correspond to otrr^(mz) = 0.115 and 0.125 respectively. 



The lattice spacing is determined from 1S-1P and 
1S-2S quarkonium level splittings which are known 
experimentally to be very insensitive to the heavy 
quark mass in the bottom to charm region. Quarko- 
nia are tiny systems, so finite volume effects should 
be small. The actions used can be systematically im- 
proved to control discretisation errors. Unquenching 
effects can be estimated by extrapolating Qp 1 linearly 
innf. A further systematic error is the effect of the sea 
quark mass on the level splittings. Sea quarks in lat- 
tice simulations are heavier than physical up or down 
quarks. Grinstein and RothsteinO estimate that the 
extrapolation to physical masses could increase the 
final value for oty^a{mz) by 0.003. 

To convert to a continuum coupling, one uses: 



7.1 a s from Quarkonia 

Tta_measured lattice coupling is ap, defined exactly 
byE 

47r /3.4 W r , s 

-lnWi,i = y a p(— )(!- [l-19 + 0.07n/]a P ). 

Wj.,1 is the single plaquette expectation value which 
can be determined accurately for n./=0 and 2 with 
varying sea quark masses. 



c£g(Q) = a<?' ) (Qe*'*)(l + la P 



MS 



+ •••)• 



The n /-dependent constant C2 is now known III in 

the quenched case, Cj"' -0 '' = 0.96. ft was previously 
set to zero. Using the quenched value even for nj = 3 
raises the result for a-^(mz), although there is clearly 
still a systematic error here. 

Fig. |^ shows recent measurements of ap, plotted 
as ap 1 against the value Q = 3.41/a at which they 
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are measuredEZI The solid curves superimposed on the 
figure show what value for «p 1 (Q) for n/ = 3 would 
convert to given values for a^(mz). The moat recent 
results from the NRQCD collaboration givel£3: 



o- 



MS 



(mz) 



0.1175(11)(13)(19) 1S-1P, 
0.1180(14)(14)(19) 1S-2S. 



The first error is a combination of statistics, determi- 
nation of the lattice spacing a and relativistic correc- 
tions. The second error comes from the extrapolation 
in the sea quark mass and the third, dominant, error 
is from the conversion to (from the difference be- 
tween using C2 = and 0.96). These combine to give 
a m {m z ) = 0.118(3). The NRQCD method currently 
relies on perturbation theory to fix the coefficients in 
the action, for which it is hard to estimate systematic 
errors. Nonperturbative renormalisation techniques 
may help here. The Fermilab/SCRI group use a dif- 
ferent effective action to calculate the bb and cc spec- 
tra in the quenched approximation and the bb spec- 
trum for nj=2. Applying tiie same analysis as above 
yields a preliminary resultEll of a^g(mz) = 0.116(3). 
Combining the NRQCD and Fermilab/SCRI results I 
quote 



W = 0.117(3). 



8 Light Quark Masses 

The masses of the light quarks, m u , m<i and m s , are 
three of the least well known standard model param- 
eters, but their values are important in a number of 
areas. The strange quark mass, for example, appears 
in the evaluation of matrix elements for the AI =1/2 
rule and the CP-violation parameter e'/e. Chiral per- 
turbation theory allows the extraction of the mass ra- 
tios from pseudoscalar meson masses. QCD sum rules 
can be applied for the masses themselves, but rely on 
detailed experimental information about the hadronic 
spectral function. Direct calculation of the masses 
from lattice QCD is thus an important challenge. 

Lattice simulations determine a bare lattice- 
spacing-dependent quark mass m(a) which can be re- 
lated to a continuum renormalised mass, such as the 



MS mass, 



,MS 



(/i). The conversion factor can be cal- 



culated using (boosted) perturbation theory. It has 
become standard to quote results for ra MS (2 GeV), us- 
ing a scale which matches the typical scale of lattice 
calculations (a -1 ~ 2-4 GeV) and for which pertur- 
bation theory should work better. The lattice mass is 
determined by evaluating pseudoscalar or vector me- 
son masses, for which the mass-squared or mass itself 
depend linearly on the quark mass respectively. 

Gupta and Bhattacharya (G&B)Oo have made 
a recent global analysis, performing a continuum ex- 




Figure 6: Strange quark mass, m^ ls (2GeV), extracted us- 
ing m0i-m-T«uenched calculations with Wilson and staggered 
quarks LJ'La Also shown are fits for the continuum extrapola- 
tion of both sets of points. 



trapolation, a — > 0, of quark masses calculated us- 
ing both Wilson and staggered fermions. The ex- 
trapolation for the strange quark mass is shown in 
Fig. ||. For Wilson quarks the leading lattice spac- 
ing dependence is 0(a) and for staggered quarks 
it is 0(a 2 ). A_Chicago-Fermilab-Hiroshima-Illinois 
(CFHI) groupOL3 compare results for Wilson quarks 
and quarks with a tadpole-improved SW action, where 
O(a) effects are reduced, first appearing as 0(a s a), 
and also perform a continuum extrapolation. 

Results for the strange quark mass from quenched 
calculations are: 



mf=(2GeV) 



f 90(20)(10)MeV G&B^_ 
\ 95(16) MeV CFHlE 



For the u and d quarks, the quantity extracted 



is the average mass, m — (m u +md)/2. Quenched 
results for this are: 

3.2(4)(3)MeV G&B@i 
3.6(6) MeV CFHlE 



m Mb (2 GeV) 



These values are already at the bottom end of the 
range predicted by other methods. Since m s appears 
quadratically in the evaluation of matrix elements for 
e'/e, a low value could have important implications 
for standard model calculations of CP violationEJ 

The Wilson data for the strange quark mass 
clearly depend on the lattice spacing. Although the 
leading a dependence is linear, it could be that a num- 
ber of effects are conspiring to produce an apparent 
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linear dependence in the current data (compare the 
case for Bk for staggered fermions, where the lead- 
ing 0(a 2 ) is satisfied, if at all, only for the data at 
the smallest lattice spacings used to date). In this 
case, the final result may well turn out to be differ- 
ent. These values rely on a perturbative matching be- 
tween lattice and continuum definitions of the quark 
masses, which could prove to be unreliable although 
the perturbative correction is not large. Nonpertur- 
bative methods have been proposed and tested, based 
on the Ward identity for the axial vector current EZl 

The results above are all in the quenched approxi- 
mation. Some calculations are also available with two 
flavours of dynamical fermions. For the same lattice 
spacing they lie systematically below the quenched re- 
sults, an effect anticipated from the different running 
of the_strong coupling in quenched and unquenched 
QCDL3 The reduction of the quark masses in full 
QCD might be very dramatic, but I believe the current 
paucity of data forbids meaningful numerical predic- 
tions. It will be extremely interesting, however, to 
follow developments in these calculations. 

9 The Lightest Scalar Glueball 

Experiments El show that there are more scalar- 
isoscalar resonances with masses below 1.8 GeV than 
can be accounted for by light qq states. This is strong 
evidence that glueballs exist and two pxime candidates 
are the / (1500) and / O (1720) mesonsEl Since glueball 
and qq states are expected to mix, the glueball content 
of these candidates remains to be determined. 

Lattice calculations provide input by calculating 
glueball masses in the quenched approximation, where 
glueballs are stable and do not mix with quark states. 
In principle, full QCD can be modelled on the lattice, 
tuning the sea quark masses between physical values, 
where the experimental meson spectrum should be re- 
produced, and large values, where one matches onto 
quenched results. In practice, glueball studies using 
dynamical quarks, where glueball-meson mixing can 
occur, are just beginning. Initial results Ell suggest 
that the masses will not change by more than about 
10%. The dynamical masses used are still quite large, 
however, and things will become more complicated 
when these masses are light enough to allow glueball 
decay. 

Two groups have continuum- 

extrapolated quenched results for the lightest scalar 
glueball: 

_ / 1550(50) MeV UKQCTjIl 
m ° ++ ~ \ 1740(71) MeV GF11E2 

Continuum extrapolations of all available data can dif- 



fer because different quantities are used to fix the lat- 
tice spacing and different ranges of a may be used. 
Results are: 

( 1600(160) MeV MichaelHIL 
m ++ = < 1707(64) MeV WeingartenS 

[ 1610(150) MeV Close & TeperBH 

Quenched lattice.jcalculations also now exist for scalar 
ss quarkoniumEZrEa These have not been extrapolated 
to zero lattice spacing, but the evidence is that the 
quenched scalar ss mass is below the quenched scalar 
glueball mass. 

The lattice results can be used as input Ear 
glueball-quarkonium mixing models. Weingarten 
has a simple model mixing the glueball with ss and 
nn = {uu + dd)/V2. The observed / (1390), /o(1500) 
and /o(1720) masses (if the latter state is confirmed as 
a scalar) can be reproduced for input glueball and ss 
masses of about 1640 MeV and 1520 MeV respectively, 
consistent with the lattice results. The /o(1720) is 
more than 75% glueball (in probability), while the 
/o(1500) is more than 75% ss. Moreover, the nn 
component of the /o(1500) has opposite sign to the 
ss component, which could help explain the observed 
suppression in the width of /o(1500) decays to KK. 

A quenched calculation has also been made of 
the coupling of the ++ glueball to two pseudoscalar 
mesonsEa This is a very delicate calculation made 
at a single value of the lattice spacing. The result, 
however, is a value of 108(29) MeV for the total two- 
body width, implying that the lightest scalar glueball 
should be easy to find experimentally. The calcula- 
tion also indicates that the two-body couplings in- 
crease with increasing pseudoscalar mass, consistent 
with observations for the /o(1720). 

The Weingarten model and arguments sug- 
gest that the / (1720) is predominantly a glue- 
ball and /o(1500) is predominantly ss quarko- 
nium. However, other models produce different 
conclusions E2k3 Various experimental tests to deter- 
mine the flavour rotitept of the isoscalar mesons have 
been proposedtSE3o'E3 Future lattice calculations, 
models and experiments should help pin down the 
glueball. 
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